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Abstract 

System-centric modeling and analysis are of key significance in planning and optimizing cellular networks. 
In this paper, we provide a mathematical analysis of performance modeling for LTE networks. The system model 
characterizes the coupling relation between the cell load factors, taking into account non-uniform traffic demand 
and interference between the cells with arbitrary network topology. Solving the model enables a network- 
wide performance evaluation in resource consumption. We develop and prove both sufficient and necessary 
conditions for the feasibility of the load-coupling system, and provide results related to computational aspects 
for numerically approaching the solution. The theoretical findings are accompanied with experimental results to 
instructively illustrate the application in optimizing LTE network configuration. 

Index Terms 

3.5G and 4G technologies, cell load coupling, network planning, optimization, system modeling 

I. Introduction 

Planning and optimization of LTE network deployment, such as base station (BS) location and 
antenna parameter configuration, necessitate modeling and algorithmic approaches for network-level 
performance evaluation. Finding the optimal network design and configuration amounts to solving 
an optimization problem of combinatorial nature. Toward this end, system modeling admitting rapid 
performance assessment in order to facilitate the selection among candidate configuration solutions, of 
which the number is typically huge, is essential. In this paper, we provide a rigorous analysis of an 
LTE system performance model that works for general network topology and explicitly accounts for 
non-uniform traffic demand. The performance model that we study is referred to as the load-coupling 
system, to emphasize the fact that the model characterizes the coupling relation between the cells in 
their load factors. For each cell, the load factor is defined as the amount of resource consumption in 
relation to that is available in the cell. The load value grows with the cell's traffic demand and the 
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amount of inter-cell interference. Intuitively, low load means that the network has more than enough 
capacity to meet the demand, whilst high load indicates poor performance in terms of congestion and 
potential service outage. In the latter case, the network design and configuration solution in question 
should be revised, by reconfiguration or adding BS infrastructure. Thus, simple means for evaluating the 
cell load for a given candidate design solution is of high importance, particularly because the evaluation 
may have to be conducted for a large number of user demand and network configuration scenarios. 

The load-coupling model for LTE networks takes the form of a non-trivial system of non-linear equa- 
tions. Calculating the solution to the model, or determining solution existence, is not straightforward. In 
this paper, we present contributions to characterizing and solving the load-coupling system model. First, 
we present a rigorous mathematical analysis of fundamental properties of the system and its solution. 
Second, we develop and prove a sufficient and necessary condition for solution existence. Third, we 
provide theoretical results that are important for numerically approaching the solution or delivering a 
bounding interval. Fourth, we instructively illustrate the application of the system model for optimizing 
LTE network configuration. 

The remainder of the paper is organized as follows. In Section |II] we review some related works. The 
system model is presented in Section Unl and its fundamental properties are discussed in Section |IVl 
In Section |Vl we present linear equation systems for the purpose of determining solution existence. In 
Section |Vll we provide the relation between solving the load-coupling system and convex optimization, 
and discuss approximate solutions. The application of the system model and our theoretical results to 
LTE network optimization is illustrated in Section IVIIl and conclusions are given in Section IVIIIi 

II. Related Works 

Planning and performance optimization in cellular networks form a very active line of research in 
wireless communications. There are many works on UMTS network planning and optimization. The 
research topics range from BS location and coverage planning JH-jSl, [|25l . Il43l . antenna parameter 
configuration ifTSl . lfT6l . Il33l . to cell load balancing ifTSl . Il34l . For UMTS, the power control mechanism 
that links together the cells in resource consumption is an important aspect in performance modeling O, 
lEl, lfT9l , Bm, [|42ll . By power control, the transmit power of each link is adjusted to meet a given signal- 
to-interference-and-noise ratio (SINR) threshold. By the SINK requirement, the power expenditure of 
one cell is a linear function in those of the other cells. As a result, the power control mechanism is 
represented by a system of linear equations, which sometimes is referred to as UMTS interference 
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coupling IfTSl . |fT6l . Interference coupling can be modeled for both downlink and uplink. For network 
planning, the interference coupling system needs to be solved many times for performance evaluation 
of different candidate network configurations and multiple or aggregate user demand snapshots. In [|26ll . 
it is shown that, for both downlink and uplink radio network planning, the dimension of the power- 
control-based system of equations can be reduced from the number of users in the system to the number 
of cells. The observation stems from system characteristics that also form the foundation of distributed 
power control mechanisms, see, e.g., |fT9ll , ||42]| . In lfT2l . the authors provide theoretical properties 
of the power-control-based system, and feasibility conditions in terms of target data rates and QoS 
requirements. Motivated by the fact that full-scale dynamic simulation is not computationally affordable 
for large networks, the authors of H4l extend the UMTS power-control system by a randomization-based 
procedure of service and rate adaptation for HSUPA network planning. 

In cellular network planning, the power-control equation system is considered under given SINR 
threshold. Thus the system solution and its existence are induced by the (candidate) network config- 
uration in question. In a more general context of wireless communications, power control is often a 
means in performance optimization, that is, the powers are optimization variables in minimizing or 
maximizing objectives representing error probability, utility, QoS, etc., that are all functions of SINR. 
There is a vast amount of theoretical analysis and algorithmic approaches for power optimization under 
various (typically non-linear) objective functions, where a gain matrix defines interference coupling 
Il371 - B0l . In [[37ll , the authors identify objective functions admitting a convex formulation of power 
optimization, and develop a distributed gradient-projection-based algorithm. Further developments in- 
clude algorithmic design utilizing Kuhn-Tucker condition ||39ll , conditional Newton iteration yielding 
quadratic convergence [|40l , and model extension to include explicit SINR-threshold constraints [[381 . 

Another line of research of power control is the characterization of the achievable performance region 
under various utility functions and interference functions. The authors of [[TTI show the strict convexity 
of the region for logarithmic functions of SINR. In [|7][, the authors characterize utility functions and 
function transformation of power, for which the resulting power optimization problem is convex. The 
investigation in ^ provides conditions under which the boundary points of the region are Pareto- 
optimal. In [[51, the authors present graph representations of power and interference, and study the 
relation between graph structure, irreducibility of the interference coupling matrix, and the convexity 
of the utility region. 
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In contrast to the power-control model, the service requirement of rate-control scheme in cellular 
networks is not SINR threshold (or a function of SINR), but the amount of data to be served over 
a given time period. Among other advantages, this approach makes it possible to capture the effect 
of scheduling without the need of explicitly modeling full details of scheduling algorithms. The rate- 
control-based approach is primarily targeting, although not limited to, non-power-controlled systems or 
systems with a fixed-rate traffic demand. The approach has been less studied, but is of a high interest 
for OFDMA-based networks. In general, the rate-control scheme exhibits non-linear relations between 
the cell-coupling elements (in our case, cell loads). The resulting model is therefore more complex 
than the power-control model for UMTS. For power control, fundamental solution characterizations 
are well-established for linear as well as more general interference functions. For the latter, see, for 
example, HOl. For rate-control-based coupling systems (see ll23l and Section HIH). a structural difference 
from power control is that, in the former, one element cannot be expressed as a sum of terms, each 
being a function denoting the impact of another element, and the coupling is not scale invariant. For 
network planning, one known approach is to consider an approximate linear function, obtained from 
system-specific adaptive modulation and coding (AMC) parameters, to represent the relation between 
date rate and SINR ll24ll . and thereby arrive at a equation system being similar to that of UMTS. 

From an engineering standpoint, UTE network optimization is becoming increasingly important. In 
lfT3l . the authors provide the fundamental principles of UTE network operation and radio resource 
allocation. Among the optimization issues, the research theme of scheduling strategies and radio resource 
management (RRM) algorithms has been extensively investigated. See, for example, |l6l, ||201 - ||22]| . [|29l - 
II3TI and the references therein. Two major aspects considered in the references are the balance between 
resource efficiency and fairness, and quality of service awareness. In IfTTl , the author gives a survey of 
tools enabling service and subscriber differentiation. For cell planning, the propagation modeling, link 
budget consideration, and performance parameters have been investigated in [|36ll . 

High-level and accurate performance modeling is of high value in planning cellular networks, as full- 
scale dynamic simulations are not affordable for large planning scenarios (e.g., Ii44|). The LTE system 
model that we analyze has been introduced by Siomina et al. [32| for studying OFDM network capacity 
region with QoS consideration. The work in ll32l does not, however, provide a general analysis of the 
model, and the major part of the study relies on a simplification assuming uniform load among cells. 
In the forthcoming sections, we present both analytical and numerical results without these limitations. 



SIOMINA AND YUAN: ANALYSIS OF CELL LOAD COUPLING FOR LTE NETWORK PLANNING AND OPTIMIZATION 5 

Recently, the authors of [|23l have presented a non-linear LTE performance model being very similar 
to the one studied in the current paper. That our performance model has been independently proposed 
by others supports the modeling approach. The work in ll23ll provides further an approximation of load 
coupling via another non-linear but simpler equation system, along with incorporating continuous user 
distribution. Our study differs from ll23l . as the focus of the current paper is a detailed investigation of 
key properties and solution characterization of the load coupling system. 

in. The System Model 

Denote hy J\f = {1, ... ,n} the set of cells in a given network design solution. Without loss of 
generality, we assume that each cell has one antenna to simplify notation. The service area is represented 
by a grid of pixels or small areas, each being characterized by uniform signal propagation conditions. 
The set of pixels is denoted by J'. The total power gain between antenna i and pixel j is denoted by 
gij. We use Ji C J io denote the serving area of cell In a network planning context, both the gain 
matrix as well as the cells' serving areas are determined by BS location and antenna configuration. 

For realistic network planning scenarios, the traffic demand is irregularly distributed. Let the user 
demand in pixel j be denoted by dj. The demand represents the amount of data to be delivered to 
the users located in pixel j within the time interval under consideration. By defining a service-specific 
index, the demand parameter and the system model can be extended to multiple types of services (see 
[|32l ). We will, however, consider one service type merely for the sake of compactness. 

We use Pi to denote the level of resource consumption in cell i. The entity is also referred to as cell 
load. In LTE systems, the cell load can be interpreted as the expected fraction of the time-frequency 
resources that are scheduled to deliver data. The network- wise load vector, p = (pi, p2, • • • , PnY , plays a 
key role in performance modeling. In particular, a well-designed network shall be able to meet the target 
demand scenarios without overloading the cells. Hence the load vector forms a natural performance 
metric in network configuration (cf. power consumption in UMTS networks). The load of a cell is a 
result of the user demands in the pixels in the cell serving area, the channel conditions, as well as the 
amount of interference. The last aspect inter-connects the elements in the load vector, as the load of 
a cell is determined by the SINRs and the resulting bit rates over the cell's serving area, and these 
values are in turn dependent on the load values of the other cells. To derive the performance model, 
we consider the SINR in pixel j G Ji defined as follows. 
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iAp) = -^^^ (1) 

2^ Pkgkjpk + cr 
keM\{i} 

In ([I]), Pi is the power spectral density per minimum resource unit in scheduling (in LTE, this 
corresponds to a pair of time-consecutive resource blocks), and cr^ is the noise power. By ([T]), the inter- 
cell interference grows by the load factor. In effect, pk can be interpreted as the probability of receiving 
interference originating from cell k on all the sub-carriers of the resource unit. Let B log2(l +7j(p)) be 
a function describing the effective bitrate per resource unit. This formula is shown to be very accurate 
for LTE downlink Il271 . Thus to serve demand dj in j, ^i^g (i!j,^.(p)) resource units are required. 

Let K denote the total number of resource units in the frequency-time domain in question, and denote 
by pij the proportion of resource consumption of cell i due to serving the users in j E Ji. By these 
definitions, we obtain the following equation. 

From (|2l), it is clear that the load of a cell is a function of the load levels of other cells. Observing 
that Pi = YjjeJi Pv ^^'^ putting the equations together lead to the following equation. 

The equation above represents the coupling relation between cells in their resource consumption. 
In vector form, we have p = f{p, g, d, K, B), where / = (/i, . . . , /i, . . . , and /», i = 1, . . . , n, 
represents the R"^^ R+ function as defined by ([3]); here, and are used to denote the single- 
and {n — 1) -dimension space of all real non-negative numbers, respectively. Since in the subsequent 
discussions there will be no ambiguity in the input parameters, we use the following compact notation 
to denote the non-linear equation system, 

P = f{p)- (4) 

From ([3]), three immediate observations follow. First, for all i = 1, . . . ,n, the load function fi is 
strictly increasing in the load of other cells. Second, for non-zero a^, this function is strictly positive 
when the load values of other cells (and thus interference) are all zeros, i.e., /(O) > 0. Third, the 
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function is continuous, and at least twice differentiable for p > 0. 

From the network performance standpoint, the capacity is sufficient to support the traffic demand, 
if equation system @) admits a load vector p with < Pi < l,i E M. In our analysis, however, we 
do not restrict p to be at most one, in order to avoid any loss of generality. In addition, even if the 
solution contains elements being greater than one, the values are of significance in network planning, 
because they carry information of the amount of shortage of resource in relation to the demand. 

Solving dH) deals with finding a fixed point (aka invariant point) of function / in M" , or determines 
that such a point does not exist. In the remainder of the paper, we use S as a general notation for the 
space of non-negative solutions (fixed points) to systems of equations or inequalities. The system in 
question is identified using subscript. Thus, Sp=fi^p) denotes the solution space of dH). Note that, for 
dH) as well as the linear equation systems to be introduced later, only non-negative solutions are of 
interest. Hence, throughout the article, a (linear or non-hnear) system is said to be feasible, if there 
exists a solution for which non-negativity holds, otherwise the system is said to be infeasible (even if 
a solution of negative values exists). The case that dH) is infeasible is denoted by S'p=/(p) = 0. 

A useful optimization formulation in our analysis is the minimization of the total cell load, subject 
to the inequality form of dl]). The formulation is given below. 

min ^ pi (5a) 

P > /(P) (5b) 
peRl (5c) 

For dS]), its solution space Sp>f(^p) is also referred to as the feasible load region. Recall that /(p) is 
strictly increasing, hence if >S'p> 7^ 0, then for any optimal solution to dS]), (l5bl) holds with equality, as 
otherwise d5al) can be improved, contradicting that the solution is optimal. In conclusion, any optimum 
of d5]) is a solution to 

We end the section by an illustration of the load-coupling system for two cells in Figure [T] The two 
cells have symmetric parameters. In the figure, the two non-linear functions are given by the solid lines. 
In the first two cases, system dH) has solutions in M^, though one of them represents a solution beyond 
the network capacity. In the last case, the system is infeasible, as the two curves will never intersect 
in the first quadrant. The straight lines with markers in the figure represent linear equations related to 
dl]). Details of these linear equations are deferred to Section |Vl 
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IV. Fundamental Properties 

In this section, we present and prove some fundamental properties of the load-coupling system (HI). 
These theoretical results are of key importance in the study of solution existence and computation. 
For compactness, we introduce additional notation to simplify (|3]) while keeping the essence of the 
equation. Define aj = bikj = -p^, and Cij = -p^- These parameters contain, respectively, the 
relation between the demand in pixel j and the resource in cell i, the inter-cell coupling in gain between 
cells k and i in pixel j, and the channel quality of cell i in relation to noise in pixel j. The load equation 
^ can then be written in the following form. 

Pi = fi{p), where fi{p) = V — — ^ -. (6) 

a,- logo 1 + \ ; — 

The first fundamental property of ([4]) is how fast the load of a cell asymptotically grows in the load 
of another cell. We formulate and prove the fact that, in the limit, the first-order partial derivative of 
the load function converges to a constant. For any two cells i, k (i ^ k), |^ is equal to 



^\n(2)^— J i - J . (7) 



Theorem 1: lim — ^ = y^ln(2)-^ 

Pfe^oo dpk ^ ttj 

Proof: Consider the component for pixel j in the sum in ([T]), and ignore the constant multiplier 
ln(2)^. Letting u = Xl/iGA^\{j} ^ihjPh + Cij, (|7]) can be written as the following expression. 



u\l + I) ln(l + i) ln(l + i) ln(l + 1)- ln(l + + ln(l + 1)« ln(l + i) 

The theorem follows then from the facts that lim„_^.oo(l + = ^ ^i^d u is linear in p^. ■ 

By Theorem [H the load of a cell increases linearly in the load of another cell in the limit, i.e., the 
function converges to a line in the high-load region. Moreover, the slope of the line is strictly positive. 

The next fundamental property is concavity. The examples in Figure [T] indicate that the load of a 
cell is a strictly concave function in the other cell's load. We show that this is generally true. 

Theorem 2: For any cell i G M, fi is strictly concave for (pj, . . . , p^+i, . . . , p„) G 
Proof: Without loss of generality, consider and its (n — 1) x (n — 1) Hessian matrix. Let 
u = J2h£Af\{n} ^nhjPh + Cnj. FoT two cells k and h, the Hessian element has the following expression. 
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Let q{u) = 2 — (2u + l)ln(l + We show that q{u) < for m > 0. This holds, for example, 
for u = 1. Next, lim^^o, ^(m) = 2 - lim™ In [(1 + + + ^)] = 0. Consider q'{u) and 
q'{u) = -2 In (l + i) + i + and ^"(m) = ^.^'l^y, < 0,Vm > 0. Therefore ^'(m) is strictly 
decreasing and limu^oo = 0. Hence, q'{u) > 0, meaning that q{u) is strictly increasing for u> 0. 
This, together with g(l) < and lim„_^oo = 0, prove that q{u) < 0, Vn > 0. By the definition 
of u, for p > 0, u > Cnj which a strictly positive number for non-zero noise power. Hence ([8]) is 
well-defined and negative for all p > 0. Next, observe that the Hessian matrix is the result of the 
following expression. 

Y ln(2)ln(l + i)g(u) 
ftj. «^ [\n\l + l){u^ + u)f 

Because of the form of ^ and that q{u) < 0,Vm > c„j > 0, the Hessian matrix is negative definite 
for any p > 0. Hence the conclusion. ■ 

From the concavity result, it follows that, for any cell i, . . . , pj_i, pj+i, . . . , p„) — pi exhibits 

a strict radially quasiconcave structure. A function is radially quasiconcave, if for a given stationary 
positive point, which is in our case a solution to /i(pi, . . . ,pj„i,pj+i, . . . ,p„) = pi, and any scalar 
in range (0, 1), the function value of the scaled point is greater than or equal to zero. If the value is 
positive, the function is strictly radially quasiconcave. 

Corollary 3: For each i G M, /i(pi, . . . , Pi-i, Pi+i, . . . , Pn) — Pi is strictly radially quasiconcave, i.e., 
if /i(pi, . . . , pi_i, pi+i, . . . , p„) = Pi, then /i(Api, . . . , Api„i, Api+i, . . . , Ap„) > Ap^ for any A G (0, 1). 

Proof: Note that /i(A(pi, . . . , pi_i, pi+i, . . . , p„)) = /i(A(pi, . . . , pi_i, p^+i, . . . , p„) + 0(1 - A). By 
Theorem H we have fi{\{pi, Pi_i, p^+i, . . . , p„)) > A/i(pi, . . . , pi_i, p^+i, . . . , p„) + (1 - A)/i(0). 
Since /i((pi, . . . , Pi-i, Pi+i, . . . , p„)) = pi and /^(O) > 0, the result follows. ■ 

In a real-life LTE network, if the capacity is sufficient to accommodate the demand, then the network 
load will be at a stable working point, which should be unique. Thus the performance model / is 
reasonable only if uniqueness holds mathematically. The following theorem states this is indeed the 
case. In the rest of the paper, the unique solution, if it exists, is denoted by p*. 

Theorem 4: If S'/(p)=p ^ 0, then it is a singleton, i.e., /(p) = p has at most one solution p* in M". 
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Proof: Suppose there are two solutions and p^, both satisfying Q, and p^ ^ p^. Let m E 
argmin-^^ nPl/pf, and A = p^/ Pm- Thus p}^^ = Ap^. Assume A < 1. Then by construction, Ap^ < p^, 
and because / is strictly increasing in the domain of M", /,„(Ap^) < fm{p^)- Also, by Lemma [3l 
^Pm < /m(Ap^), and thus /„,(p^) > Ap^. Note that /„,(p^) = = Ap^ gives an contradiction. 
Therefore A > 1 . Considering scaling down p^ with A instead, and applying the same line of argument, 
a similar contradiction is obtained. Hence the conclusion. ■ 

V. Determining Solution Existence and Lower Bounding 

Having proven solution uniqueness, we examine the existence of p*, that is, whether or not (Hj) has a 
fixed point. There are a number of theorems characterizing the existence of a fixed point (e.g., Brouwer's 
fixed-point theorem in topology). However, these results do not apply to @ because, in general, the 
output of function / is not confined to a compact set in R" . In this section, we use a linear equation 
system for analyzing solution existence. To this end, we first present and prove some basic properties 
of the optimization formulation ([5]). 

Theorem 5: Assume Sp>f[p) ^ 0, i.e., there exists p > f{p) > 0, then ^ has an optimal solution. 
Proof: Consider the optimization problem min^,.g^pj,p G S, where S = Sp>f{p) fl {p < p}. 
By the assumption in the theorem, S* 7^ 0. From the definition of S, it is clear that any point being 
arbitrarily close to S (i.e., boundary point) is in the set, thus S is closed. In addition, S is bounded 
since S C n {p < p}. Hence S is compact, and the result follows from Weierstrass theorem in 
optimization. ■ 

Corollary 6: If there exists p > /(p) > 0, then S'p=/(p) 7^ 0. 
Proof: Follows immediately from Theorem [5] and the previously made observation that any optimal 
solution to © satisfies (l5bl) with equality. ■ 

To further characterize solution existence, we define the following type of linear equation systems, 

p = h{p) = Hip-p) + f{p), (10) 

where h = {hi, ... , hn) is a vector of linear functions, and each of them is defined in W^T^ — t- M+, 
p is a vector in M" with given values, and /(p) is a vector- function with elements defined by In 
(fTOl) . i? is an n X n matrix where the diagonal elements, Ha^i = 1, . . . ,n are zeros, and the other 
elements Hii^,i 7^ k, are strictly positive. Note that if H is the Jacobian of function / evaluated at 
point p, (fTOl) is a linearization of the non-linear equation system (H]) where the right-hand side of 
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(fTOl ) represents the tangent hyperplane to function /(p) at p. Such linear approximations are further 
discussed in Section |VIl 

Observing the fact that the partial derivative ^ asymptotically approaches a constant, as formulated 
in Theorem (TJ we consider linear approximation of / by means of the linear function having the limit 
values of the partial derivatives as the matrix elements in H, and passing through the point defined by 
the load function values with zero load. Define the case of h where = ln(2) ^ikj/o^j for 

k ^ i, and p = 0. For this linear approximation, there are similarities between the elements of H and 
the UMTS interference-coupling matrix (see, e.g., ifTSl . Il28ll ) in that both capture the relation between 
gain factors of the serving and interfering cells; however, the target QoS in the interference-coupling 
matrix is link quality, whilst in H it is given by the amount of user traffic demand. 

If Sp^f^O(^p^ ^ 0, the solution, denoted by p^, is clearly unique. The lemma below states that the 
linear function provides an under-estimation of the true load function /, thus p^, if exists, gives a 
lower bound on the solution to the non-linear system (Hj). 

Lemma 7: hP{p) < f{p) for any p > 0. 

Proof: We prove the validity of the result for an arbitrary cell i, that is, h'^Xp) < fi{p), p > 0. 
Because both /i°(p) and /i(p) are formed by a sum over j G Ji, it is sufficient to establish the inequality 
for any j E J.^. Let u = YjkeM\{i} ^ikjPk + Cij- The proof boils down to showing the following inequality. 

- - C) .n(2) . ,n(2, - - c,,)) > j-J^ ^ ^ 

Note that u > Cij by definition. The inequality holds as equality for u = Cij. It is then sufficient to 
prove that _^ i — (u — Cij) is increasing for u > Cy. Taking the derivative and doing some simple 
manipulations, one can conclude that the derivative is non-negative corresponds to the inequality below. 

= + 1) W + < 1, u > Cij (11) 

One can show easily that lim„^o+ q{u) = 0, hence (fTTI) is satisfied for some u < Cij. Moreover, 
\imu^^ q{u) = lim„_,ooln(l + ^)"ln(l + i)" + ln(l + i)"ln(l + ^) = 1. Hence it suffices to prove 
that q'{u) = (2m + 1) In^ (l + i) - 2 In (l + i) > 0, m > 0. Using the fact that ln(l + ^) > for all 
u > 0, the non-negativity of q'{u) for u >0 becomes equivalent to that the second numerator in ^ is 
negative, which is proven in the proof of Theorem [21 and the result follows. ■ 
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From Lemma Ul one can expect that the load-coupling system @i has a solution, only if a solution 
exists to p = hP{p). The following theorem formalizes this necessary condition, and establishes the 
result that bounds p* from below. 

Theorem 8: If 5p=j(p) ^ 0, then Sp^^o^^^^ ^ and pO < p*. 
Proof: Consider the following linear programming (LP) formulation. 

min ^^Pi (12a) 

P > h%p) (12b) 
peRl (12c) 

Similar to the result in Theorem [51 it can be easily proven that (fT2l) has an optimal solution if there 
exists any p > satisfying (I12bl) . In addition, it is clear that any optimum is in S'^^^o^p), and S'^^^o^p-, 
is either empty or a singleton. Consider p*. By Lemma|7l hP{p*) < f{p*) = p* ■ Hence p* is a feasible 
solution to (fT2l) . It follows then Sp^y^o^^^-^ ^ 0. Furthermore, (fT2l) obviously remains feasible with the 
additional constraint p < p* . Since the LP optimum is unique and equal to p^, p^ < p*. ■ 

By Theorem [8l the linear system p = hP{p) is potentially useful for detecting infeasibility. If the 
linear system is infeasible, then it is not meaningful to attempt to solve (H)). In addition, if feasibility 
holds for p = hP{p), the solution provides a lower bound to the true load values. Thus having p^ 
close to one indicates an overloaded network, and its corresponding configuration can be discarded 
from further consideration in network planning, without the need of solving the non-linear system 

In Figured! the lines with markers represent the linear function . In the first two cases, p* exists, 
and solving the linear system leads to a lower bound p^ (i.e., the intersection point of the lines) of p*. 
In the last case, the linear system has no solution, and consequently Sp=f(p) = 0. 

Thus far, it has become clear that p = h^{p) provides an optimistic view of the cell load. We are able 
to prove a slightly unexpected but much stronger result. The linear equations p = hP{p), in fact, give 
an exact characterization of solution existence of the load-coupling system. Namely, that p = hP{p) 
has a solution is not only a necessary, but also a sufficient condition for the feasibility of dU. 

The intuition of the sufficiency result is as follows. Consider Figures |l(a)Hl(b)| for which the linear 



equation system has solution. Suppose the slopes of the lines are increased slightly. Intuitively, if the 
increase is sufficiently small, the new linear system will remain feasible. Also, the figure gives the hint 
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that the modified linear function will eventually go above the non-linear load function for large load, 
indicating Sp=f(^p) ^ 0. To rigorously prove the result, we define the linear equation system p = h^{p), 
obtained by increasing the slope coefficients of by a positive constant e. That is, h'^ denotes the 
case of ( flOl ) where Hn. = ln(2) Xljeji hkj/O'j + e, p = 0. 

Lemma 9: If Sp^y^o^p-^ ^ 0, i.e., exists, then there exists e > such that Sp=h^p) 7^ 0- 

Proof: First, note that Sp^y^Of^p-^ has a non-empty interior. In particular, it is easily verified that 
Ap^ is an interior point for any A > 1. Denote by p such a point, that is, pi > h'^{p),i E J\f. 
Letting = ^ p. = h^^(^p) + e,^^ pfc,z G J\f. Next, set e = miUi^^ei. Then pi > 

^i'(p) + '^J2k£j\f\{i} Pki i ^ Thus for this value of e, p G Sp>h^p), and the result follows. ■ 

Lemma 10: Consider any p > and any e > 0. Denote by A a positive number. For any i G M, 
lim [h\{\p)-U{\p)] = oo. 

A— >oo 

Proof: Consider the definitions of hl{\p) and fi{\p). After some straightforward re- writing and 
ignoring the constant term /(O) in hl{\p), the difference between the two functions has the following 
form. 



^ ln(2) 



bikjPk^ 



+ 1 E E^Ah (13) 

.ke^^\{i}j€J^ ^ 



Let g(A) denote the expression in the square brackets of (fT3l) . By repeatedly using I'Hopital's rule, 
one can show that limA^oo Q'(A) = — | — Qj, which is a constant. Observing that the last term in (fT3l) 
grows linearly in A, the lemma follows. ■ 

Theorem 11: If S'p^^O(p) 7^ 0, then Sp=f(p) ^ 0. 
Proof: By Lemma [9l there exists e > and p^ satisfying p^ = h^{pl^). It is easily verified 
that Ap^ > h,'(Ap^),A > L Using Lemma [M there exists A such that Ap^ > /i'(Ap^) > fi^Ph)- 
Therefore Sp>fi^p) ^ 0, and the result follows from Corollary [6l ■ 

Theorems [8] and [m together provide a complete answer to the solution existence of LTE load coupling, 
that is, whether or not the system has a fixed point in is equivalent to the feasibility of the linear 
equation system p = hP{p). Clearly, given an LTE network design, this feasibility check should be 
performed first, before determining the load values. Furthermore, from Theorem [8l violating p^ < 1 is 
a simple indication of that p* is beyond the network capacity. For a two-cell example, the solution to 
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the linear system p = hP{p) is 

i — _M21-ni2 
i — /i2l-ni2 

With (fT4l) . a feasible solution exists when \ — H2iHi2 > 0, i.e., H12 = forms the (open) boundary 
of the feasibility region in the two coefficients. Note that H12 and H21 are linear in the traffic demands 
to be satisfied in cell 1 and cell 2, respectively. The derived relation representing the resource sharing 
trade-off for the two neighbor cells in this example is well in line with the commonly known radio 
resource sharing and capacity region concepts. 

VI. Convex Optimization and Upper Bounding 

Provided that 5'p=^(p) 7^ 0, a solution algorithm needs to be applied to find p*. Solving p = f{p) is 
equivalent to finding the (unique) root of the n-dimensional function p — f{p). Thus one approach is 
to use the Newton-Raphson method. In this section, we show that approaching p* can alternatively be 
viewed as the convex optimization problem formulated below. 

max ^^Pi (15a) 

ieJV 

p-fip)<0 (15b) 

p G M+ (15c) 

Corollary 12: Formulation (flSl) is a convex optimization problem, and if Sp=f(^p) 7^ 0, then p* is the 
unique optimum to (fTSl ). 

Proof: Because /(p) is concave (Theorem [21), p— /(p) is convex in R". Thus p— /(p) < is a 
convex set. The proof is complete by observing that, similar to Q, optimum to (fT?)) must satisfy (|15bl) 
with equality, and p* is the unique solution to p = /(p). ■ 

Following Corollary [121 any convex optimization solver can be used to approach p*. In network 
planning, one will need to solve © repeatedly to evaluate many candidate BS location and antenna 
configurations. Typically, the performance evaluation does not have to be exact in order to relate the 
quality of a candidate solution to that of another. Utilizing the structure of /, we can numerically obtain 
upper bounds to p* via linear equations. Consider any p G R" . Using the partial derivatives of / at p. 
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and the point f{p), we obtain an upper approximation of / due to concavity. Formally, denote by h 
the linear function of (flOl) where p = p and Hik, defined by ©, takes the following value, 

= |^(P) = E ln(2)^ ^ 

dpk f-i a 



The positive-valued solution to the linear system p = h{p), if it exists, is denoted by p^. As 
established below, h and yield upper estimations of / and p*, respectively. 

Corollary 13: h{p) > f{p), p > 0. 

Proof: Follows immediately from the concavity of / and the definition of h. ■ 

Theorem 14: If Sp^^p) ^ 0, then Ph> p*. 
Proof: Consider the linear programming (LP) formulation max{^^g^pj : p < h{p),p G R"}. 
Similar to (fT5l) . it is easily realized that the LP formulation, if feasible, has a unique optimum satisfying 
p = h(p). Hence the unique optimum is p^. Moreover, p* is, by Corollary [131 a feasible point to the 
LP, and hence the LP remains feasible after including p > p*, and the result follows. ■ 

The process of solving the load-coupling system, e.g., an interior point method for (fT5l) . will typically 
generate a sequence of iterations approaching p* from below. By Theorem [l4l the iterations can be 
used to compute upper bounds, thus yielding an interval confining p*. In order to speed up the process 
of network optimization, performance evaluation of a candidate planning solution can use a threshold 
of the maximum size of the interval, instead of computing the exact solution of the load vector. 

Computing an upper bound pj^, involves solving a system of n linear equations. The same amount 
of computation applies to the feasibility check and computing lower bounding p^ in Section |Vl It 
is straightforward to see that calculating the coefficients is of complexity C(n^). Thus the overall 
complexity lies in the matrix inversion operation that runs in the time range 0{n'^'^™) and 0{n^), 
where the former is attainable only asymptotically by the Coppersmith-Winograd algorithm. 

VII. Numerical Results 

In this section, we numerically investigate the theoretical findings in the previous sections. An 
illustrative simulation study has been conducted for a three-site 3GPP LTE network with an inter- 
site distance of 500 m, adopting a wrap-around technique. The simulated system operates at 2 GHz 
with 10 MHz bandwidth. Each site is equipped with a three-sector downtilted directional antenna with 
14 dBi antenna gain. The propagation environment and user distribution follow the 3GPP specification 
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in |[T1, assuming propagation model 1 (Okumura-Hata, urban, 8 dB standard deviation shadow fading) 
and user generation scenario 4b with one hotspot of 40 m radius per macro cell area. Note that, as 
for any system model, the complete assessment of the model validity would also include validating of 
numerical results against results from real deployments, which is beyond the scope of the current paper. 

The network layout we have used is illustrated in Figure [2l Two layers of users are generated, with 
30 users per macro cell area in total, out of which 2/3 (the dot markers) is in a randomly placed hotspot, 
and 1/3 (the x-markers) are distributed randomly and uniformly over the area. Each user equipment has 
an omni-directional antenna with dBi antenna gain. The traffic demand corresponds to 400 kbit/s for 
all users within a duration of one second in the time domain. 

Two network configurations are illustrated in Figure [2l with the only difference being the antenna 
direction of cell 1, which impacts the sets of users served by the cell and its neighbors. Intuitively, 
configuration two is inferior, since it results in that the hotspot users in cell I's original coverage area 



(see Figure |2(a)| ) is to be served by cell 8 and/or cell 9, although these users are relatively far away 
from the two cells. The likely impact is poorer link quality for the users in the handed-over hotspot as 
well as increased number of users to be served by the neighbors of cell 1. These effects are expected 
to be seen in the load of the neighbor cells. 

First, for both configurations, the existence of system solutions has been verified by finding p° to 
the corresponding linear system, as described in Section |Vl Next, the non-linear coupling system dH) is 
solved using the non-linear optimization toolbox of MATLAB. Both p^^ and p* are shown in Figure [3l 
As expected, the load of cells 8 and 9 increases for the second configuration. At the same time, the load 
of cell 1 does not decrease either, even though it serves fewer users under the second configuration. 
This is due to a joint effect of several factors. Firstly, as can be seen from Figure |2(b)[ users served 
by cell 1 are likely to experience high interference from cell 6 and vice versa. From Figure |3(b)[ we 
observe that the load of cell 6 has also slightly increased. Secondly, the increased load in cells 8 and 9 
implies more frequent transmissions in these cells and thus higher probability of interference to other 
cells; this in turn increases the load of the other cells, which can also be clearly seen in Figure [3l We 
further note that the solution p* and thus the second configuration are not feasible from the capacity 
point of view. 

In Figure [3j we also illustrate the load solutions to the linear systems described in Section |Vll i.e., 
the upper bound pj^ assuming p = p^. Table U provides further details of the quality of both the lower 
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and upper bounds obtained for all 9 cells in Configuration 1 for which the solution is illustrated in 
Figure |2(a)[ 

The tight upper bound indicates the efficiency of the linear approximation described in Section |VIl 
In average, the estimation is deviates only a few percent from the true load value. For the values 
are significantly lower than p*, as p = h^{p) represents a very optimistic view of load coupling. The 
observation sheds further light on the importance of fundamental characterization that the two systems 
are completely equivalent in solution existence, despite the large difference in numerical values of p° 
and p* . Improving the lower bounds, although being beyond the scope of the current paper, is an 
interesting topic for future investigation. For example, in [[351 , the model discussed in this paper is 
applied for load balancing, for which very tight lower and upper bounds are obtained using few fixed- 
point iterations. It should be further noted that although the results have been presented for downlink, 
the model and the theoretical findings can be also applied to uplink. 

For the two network configurations, in Figure |4] we illustrate the behavior of the cell-load coupling 
with respect to demand, which is successively scaled up uniformly over the service area. Figure |4(a)| 
and Figure [4(b)] show, respectively, the results for the nonlinear load-coupling system dH), and the linear 
equation system p = hP{p) that provides lower estimation and characterizes feasibility of dH). The two 
configurations are distinguished by using respectively solid and dotted curves for the load solutions of 



nine cells. In Figure |4(a)[ the thicker curves represent the load of cell 8. For the linear system, only 
the maximum value among the cells is shown in Figure |4(b)| for the sake of clarity. 

From Figure [4(a)| it is apparent that the solution values of dl]) grow rapidly in the high-demand region. 



and becomes infeasible beyond some point. The feasibility boundaries for the two configurations are 
shown by the vertical lines. Configuration one is clearly superior, as its load values, shown by the solid 
curves, are below those of configuration two, and the feasibility boundary is considerably higher. For 
configuration two, cell 8 has the highest load (the dotted thick line). Using configuration one, some 
of the users are served by cell 1 instead (see Figure O, leading to lower load in cell 8 (the solid 
thick line). Note that, for both configurations, when getting somewhat close to the infeasible region, 
the solver gives solutions containing some zero elements, indicating the solution is invalid (the system 
becomes unstable), before all the values abruptly drop to zeros showing infeasibility. When this occurs, 
the distance to the feasibility boundary is, in fact, significant. The behavior shows the importance of 
our analysis of characterizing feasibility exactly by the linear system p = hP{p). 
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In Figure |4(b)[ the linear system gives values growing consistently in demand. The system gives the 
feasibility boundary point of the non-linear load coupling equations, when the determinant of / — H, 
where / is the identity matrix and H is the matrix defined for h^, equals zero. After passing the 
boundary point, the linear system returns negative (infeasible) solutions. Hence the numerical results 
verify that the significance of the linear system to identifying solution existence. Moreover, using the 
linear system, one is able to conclude, as shown in Figure H that configuration one is clearly superior 
to configuration two. 

VIII. Conclusions 

We have provided a theoretical analysis of the LTE load coupling system originally presented in 
||32l and have derived its fundamental properties, including concavity, behavior in limit, and solution 
uniqueness. We have also formulated the necessary and sufficient condition for solution existence. The 
analysis leads to a simple means for determining feasibility. In addition, we have presented two linear 
approximations. The analysis has been supported by theoretical proofs and numerical experiments and 
can serve as a fundamental basis for developing radio network planning and optimization strategies for 
LTE. Furthermore, the presented linearizations and the bounding-based optimization can potentially be 
used for more general convex optimization problems with similar properties. 
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(a) Feasible solution within network ca- (b) Feasible solution beyond network ca- (c) Infeasible system: Sp=/(p) = ' 
pacity: < p < 1. pacity: 3{i G Af)\pi > 1. 

Fig. 1. An illustration of the load-coupling system of two cells. 
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(a) Configuration 1. 
2. Network configurations for numerical studies. 
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(b) Configuration 2. 



TABLE I 

Lower and upper bounds quality for Configuration 1 



Bound 


Cell 1 


Cell 2 


Cell 3 


Cell 4 


Cell 5 


Cell 6 


Cell 7 


Cell 8 


Cell 9 


Upper bound, ■ 100% 
Lower bound, -100% 


5.19 
48.05 


5.55 
48.15 


6.66 
44.44 


4.61 
43.08 


4.93 
43.21 


4.94 
50.62 


3.33 
51.66 


5.20 
46.82 


4.17 
50.00 
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(a) Configuration 1: Feasible solution within network capacity. (b) Configuration 2: Feasible solution beyond network capacity. 
Fig. 3. Load solutions for the two example network configurations. 
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(a) Growth of cell load with respect to demand for the load (b) Growth of cell load with respect to demand for the linear 
coupling system Q. system p = h}\p). 



Fig. 4. System solution with respect to traffic demand. 



